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Abstract 



The nonlinear magneto-optical Kerr-effect (NOLIMOKE) is determined for 
freestanding Fe monolayers with several in-plane structures from first princi- 
ples. Based on the theory of nonlinear magneto-optics by Hiibner and Benne- 
mann [Phys. Rev. B, 40, 5973 (1989)] we calculate the nonlinear susceptibil- 
ities of the monolayers using the ah initio FLAPW-method WIEN95 with the 
additional implementation of spin-orbit coupling and the calculation of the 
dipole transition matrix elements appropriate for freestanding monolayers. 
We present results for the spectral dependence of the nonlinear susceptibility 
tensor elements and the resulting intensities and Kerr angles. Special empha- 
size is put on the effects of structural changes such as the variation of the 
lattice constant and different surface orientations. The influence of spin-orbit 
coupling on the tensor elements for different magnetization directions is pre- 
sented as well as the azimuthal dependence of the intensities generated by 
several low index surfaces, showing the pronounced sensitivity of second har- 
monic generation to lateral structural changes as well as magnetic properties 
even in the monolayer range. 
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I. INTRODUCTION 



The nonlinear magneto-optical Kerr-effect (NOLIMOKE) combines the surface and in- 
terface sensitivity of second harmonic generation (SHG) with the sensitivity on magnetic 
ordering, which is much stronger here than in linear optics The existence and the de- 
tectability of the nonlinear magneto-optical effect was shown independently in 1989 by Pan, 
Wei and Shen || and Hiibner and Bennemann j|] . Whereas the former work was based on a 

(2) 

group theoretical classification of the nonlinear susceptibility tensor xhU t ne latter used an 
electronic structure calculation to predict the nonlinear Kerr-effect. During the last years 
significant improvements have been made. The group theoretical analysis was extended to 
antiferromagnetic surfaces ||. Also several extensions of the electronic structure calculations 
have been performed, such as the calculation of the nonlinear magneto-optical response of Ni 
within a tight binding scheme || and of Fe layers and films using a FP-LMTO method . In 
these papers spin-orbit coupling (SOC) was treated perturbatively and only recently || SOC 
was included from first principles to obtain the nonlinear susceptibility of Fe monolayers for 
several in-plane structures. Nevertheless, no complete theory was available combining both 
the group theoretical analysis and first principles electronic structure calculations, since up 
to now the dipole transition matrix elements have not been determined from first principles. 

In this paper we present a calculation of the nonlinear magneto optics for freestanding 
Fe monolayers, which obtains all fundamental aspects of magneto-optics, the electronic 
structure including spin-polarization, SOC, and the dipole transition matrix elements from 
first principles. The symmetry aspects are completely determined by the symmetries of the 
wavef unctions. 



II. THEORY 

Our calculations are based on the formula for the nonlinear susceptibility tensor 
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derived in Therein q is the momentum of the incident light, a is a damping constant 
simulating a finite life time and chosen as 0.2 eV, f(E^i) denotes the Fermi-distribution 
and m.j is the i-th direction cosine of q with respect to the plane of incidence. The energy 
bands E^i and the wavefunctions |k, I) are determined using the ab initio FP-LAPW method 
WIEN95 0. SOC is included as described in p. Within this treatment the wavefunction 
reveal the correct symmetry of the system including SOC. Thus, by determining the dipole 
transition matrix elements (k, Z|i|k', I') from the wavefunctions, we obtain the correct form of 
the susceptibility tensor with respect to symmetry as derived in |§ from first principles. The 
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numerical differences between group theoretically allowed and forbidden tensor elements is 



about 12 orders of magnitude. The precise treatment will be described elsewhere [JIO . 

Since, at this stage of the investigations, we are mainly interested in the effects of lateral 
structural changes on the magneto-optical spectra, we calculate only freestanding monolayers 
neglecting the substrate. The absence of a substrate enables us to choose arbitrary in- 
plane structures and reduces the required computer time to a reasonable value. Of course 
the substrate cannot be neglected for a quantitative determination of the magneto-optical 
response and will be included in further studies. 

Within the electric-dipole approximation the presence of inversion symmetry breaking 
is necessary for the occurrence of SHG. Thus, in the case of freestanding monolayers, we 
have to introduce an artificial symmetry breaking perpendicular to the surface. This is 
done by integrating the dipole transition matrix elements in the upper half of the unit cell 
only, which effectively sets the wavefunctions in the lower half volume equal to a constant. 
Roughly speaking, this simulates a substrate with a constant charge density. In doing so we 
also solve the general problem of treating semi infinite systems like film or layer structures 
with a k-space method like WIEN95, which uses threedimensional translational invariance. 

The SHG intensities are then calculated from the nonlinear susceptibilities using the 
expression for the second harmonic generated field given in [IT] . 



III. RESULTS 

Fig. |l] shows the spectral dependence of the imaginary part of the four independent 
tensor elements of the nonlinear susceptibility xl% obtained from a Fe(OOl) monolayer with 
out of plane magnetization (M || z) derived from a fee lattice. The in-plane lattice constant 
is varied from 2.4 A to 2.76 A (Cu has a=2.56 A), simulating different substrates . For this 
system the nonlinear susceptibility tensor has the following form: 
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There are four independent tensor element ZCCX ZZZ - XXZ + and xyz . The first three do 
not change their sign under magnetization reversal ( + superscript) and thus will be called 
[T2| , whereas the xyz-e\emeiat changes its sign (~ superscript) and thus are 
The tensor elements show a clear dependence on the lattice constant, 

Zeros in the 



"nonmagnetic' : 



termed "magnetic". 

as was already found in calculations with constant matrix elements 
spectra are shifted to larger energies for smaller lattice constants, maxima have a larger 
magnitude for larger lattice constants. Of course, the spectral dependence becomes much 
more complicated with the inclusion of the matrix elements, since the latter affect both the 
width and absolute values of the resonances. The nonlinear Kerr angle $k p in the polar 
configuration with p input polarization is given by 
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where f s ,c are Fresnel coefficients of the incident field, A PiS and F c ^ s are transmission and 
Fresnel coefficients of the generated field and N is the refractive index of the substrate ma- 
terial, is the ellipticity and an approximation for small angles is used. From Fig. [I] it is 
clear that the contribution of the \f) z and Xxlz elements can be neglected because of their 
magnitude. Thus the nonlinear Kerr angle is directly proportional to the ratio of the Xxyz 
and the x^xx elements weighted by the transmission and Fresnel coefficients (in the numerical 
calculation, however all tensorelements are included). Though the spectral dependence of 
latter is rather weak in this energy range compared to the tensor elements, a clear relation 
between the spectral dependence of the Kerr angle and the tensor elements is not obvious. 
Nevertheless both the spectral dependence of the Kerr angles and the relevant tensor ele- 
ments for the different lattice constants (see Fig. |T|) show some common features. Maxima 
around 1 eV are followed by a crossing of the curves between 1 and 1.5 eV. The spectral 
dependence of the linear Kerr angles in the upper part of Fig. ^| shows a similar feature at 
twice the energy but exhibits much more structure for higher photon energies. More detailed 
investigations of the involved complex quantities show that the spectral dependence of the 
nonlinear Kerr angles is mainly ruled by the magnetic tensor elements. Fig. ^] shows that in 
our calculations the linear and nonlinear Kerr angles are of the same order of magnitude, 
which does not agree with the experimental findings jL3|. This is due to the specific form 



of the artificial inversion symmetry breaking performed here, which also underestimates the 
size of the zzz tensor element. Note, the present study focuses on the structural effects on 
the nonlinear Kerr spectra. This problem only occurs in the case of freestanding monolayers 
and can be overcome by including substrate effects and using unit cells for the electronic 



structure calculations in which inversion symmetry is explicitly broken flO|| . Furthermore 
the size of the tensor elements in Fig. [I] is too large compared to previous estimates [H . 
This can be a problem of the neglected screening effects at the surface. 

Fig. |3] shows the differences in the spectral dependence of the nonmagnetic tensor ele- 
ments of the Fe(OOl) monolayer with a lattice constant a=2.56 A for in-plane (M || y) and 
out of plane (M || z) magnetization. Since the two dimensional unit cell of the (001) surface 
is a square lattice, the x and y direction are equivalent for out-of-plane magnetization, and 
thus x^zlx = X z 2 yy For in-plane magnetization the x and y directions are no longer equivalent 
and thus x^Jx X z 2 yy The effect of this symmetry breaking on the nonmagnetic tensor ele- 
ments is shown in Fig. |3], yielding only a slight dependence on the magnetization direction. 
This is expected from the fact that nonmagnetic tensor elements are only of second order 
inSOCp|. 

One advantage of SHG compared to linear optics is the enhanced symmetry resolution 
due to the third rank tensor x[% involved in SHG compared to the second rank tensor xif 
in linear optics. The in-plane symmetry resolution can be monitored by calculating the 
azimuthal dependence of the intensities. This is displayed in Fig. |], where the azimuthal 
dependence of the intensities with p in andp out polarization is plotted for the Fe(001), Fe(110), 
and Fe(lll) monolayer with in-plane and out-of-plane magnetization. SHG can resolve up 
to three-fold rotational symmetries. Thus, for M || z, the SHG signal of the square lattice 
Fe(001) which has four-fold symmetry in this configuration show no azimuthal dependence. 
Also the Fe(lll) surface with M || z shows no azimuthal dependence since, in the absence of a 
substrate, this monolayer has a six-fold symmetry. The curve for the Fe(110) surface reveals 
its twofold symmetry. The dashed line indicates the intensities for inverted magnetization 
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direction, showing the effect of broken symmetries due to the magnetism and SOC. For an 
in-plane magnetization this effect is increased. The magnetization shifts the SHG intensities 
in the direction perpendicular to the magnetization. Also the well-known effect of a large 
magnetic contrast in the intensities in the transversal configuration is reproduced by our 
ab initio theory. Finally, it is clearly visible that the in-plane magnetization destroys the 
circular dependence of the intensities in the case of the Fe(OOl) and Fe(lll) monolayers. 

IV. SUMMARY AND DISCUSSION 

In summary we showed results on the nonlinear magneto optical effect of Fe monolayers 
as obtained from first principles calculations. The dipole transition matrix elements were 
computed from the wavef unctions. The inversion symmetry has been artificially broken 
to obtain SHG from the freestanding monolayers. The symmetries of the system including 
SOC are completely reflected by the wavef unctions yielding the correct form of the nonlinear 
susceptibility tensor known from group theoretical classifications ||. The tensor elements 
as well as the resulting Kerr angle show a clear dependence on the lattice constant in 
agreement with previous results using constant matrix elements. The spectral dependence 
of the nonlinear Kerr angle shows up as a complex superposition of the involved tensor 
elements and transmission and Fresnel coefficients. The precise inclusion of SOC reveals a 
weak dependence of the nonmagnetic tensor elements on the magnetization direction. The 
high sensitivity of SHG to the in-plane structure was shown for the Fe(OOl), Fe(llO), and 
Fe(lll) monolayers, showing also the important influence of the magnetic tensor elements. 

We acknowledge financial support by TMR Network NOMOKE contract No. FMRX- 
CT96-0015. 
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FIGURES 



FIG. 1. Spectral dependence of the imaginary part of the nonlinear susceptibility tensor ele- 
ments of a Fe(OOl) monolayer for different lattice constants a given in A. The magnetization is 

out-of-plane (M || z). The lattice constant is varied between 2.4 and 2.76 A. 

FIG. 2. Spectral dependence of the linear and nonlinear Kerr angle of a Fe(OOl) monolayer 

with out-of-plane magnetization. The lattice constant is varied between 2.4 and 2.76 A. 

FIG. 3. Spectral dependence of the nonmagnetic tensor elements of a Fe(OOl) monolayer with 

out-of-plane and in-plane magnetization. The weak dependence reveals the SOC dependence of 

second order. 

FIG. 4. Azimuthal dependence of the nonlinear intensities with p polarization of the incident 
and generated light for the Fe(OOl), Fe(llO) and Fe(lll) monolayer with out-of-plane and in-plane 
magnetization. 
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